Chiral superconductivity from repulsive interactions in doped graphene 
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We identify graphene as a system where chiral superconductivity can be reahzed. 
Chiral superconductivity involves a pairing gap that winds in phase around the Fermi 
surface, breaking time reversal symmetry. We consider a unique situation arising in 
graphene at a specific level of doping, w^here the density of states is singular, strongly 
enhancing the critical temperature Tc. At this doping level, the Fermi surface is nested, 
allowing superconductivity to emerge from repulsive electron-electron interactions. We 
show using a renormalization group method that superconductivity dominates over all 
competing orders for any choice of w^eak repulsive interactions. Superconductivity 
develops in a doubly degenerate, spin singlet channel, and a mean field calculation 
indicates that the superconductivity is of a chiral d + id type. We therefore predict 
that doped graphene can provide experimental realization of spin-singlet chiral super- 
conductivity. 
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Chiral superconductors feature pairing gaps that wind 
in phase around the Fermi surface (FS) by multiples of 
27r, breaking time-reversal symmetry (TRS) and parity 
and exhibiting a wealth of fascinating properties [l|-l3|- 
The search for experimental realizations of chiral su- 
perconductivity greatly intensified in the last few years 
with the advent of topological superconductivity [^ 
|6|. Here we show that chiral superconductivity with 
a d^2_y2 ± idxy {d + id) gap structure can be realized 
in graphene monolayer, a system of choice of modern 
nanoscicnce [71, Isj . We demonstrate that when graphene 
is doped to the vicinity of a Van Hove singularity in the 
density of states (DOS), repulsive electron-electron inter- 
actions induce d-wave superconductivity. Our renormal- 
ization group analysis indicates that superconductivity 
dominates over competing density wave orders, and also 
indicates that interactions select the chiral d -\- id state 
over TRS-preseving d-wave states. The nontrivial topol- 
ogy of the d + id state [l| manifests itself in exception- 
ally rich phenomenology, including a quantized spin and 
thermal Hall conductance M , and a quantized boundary 
current in magnetic field [lO[. 

The search for chiral superconductivity has a long 
history. Spin-triplet p-wavc chiral superconducti vity 
(jjx ± ipy state) has likely been found in Sr2RuOi 11| . 
which represents a solid state analog of superfluid '^He 
[l|, but the spin-singlet d-\-id state has not yet been ob- 
served experimentally. Such a state was once proposed 
as a candidate state for high Tc cuprate superconduc- 
tors [9|, |lO[, but later gave way to a more-conventional 
TRS-preserving d— wave state. The key difficulty in re- 
alizing a d-\- id state is that the interactions that favor a 
d-wave state usually have strong momentum dependence 
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How can superconductivity be induced in graphene? 
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FIG. 1: Chiral superconductivity arises when graphene is 
doped to the Van Hove singularity at the saddle point (M 
points of the Brillouin zone), a) d-\-id pairing exhibiting phase 
winding around the hexagonal FS, which breaks TRS and par- 
ity {9 = 27r/3). b) Conduction band for monolayer graphene 
[31. At 5/8 filling of the tt band, the FS is hexagonal, and 
the DOS is logarithmically divergent (c) at three inequivalent 
saddle points of the dispersion Mi (i=l,2,3). Their location is 
given by ±ei, where 2ei is a reciprocal lattice vector. The sin- 
gular DOS strongly enhances the efi^ect of interactions, driving 
the system into a chiral superconducting state (a). Since the 
FS is nested, superconductivity competes with density wave 
instabilities, and a full RG treatment is required to establish 
the dominance of superconductivity. 



Existing proposals for superconductivity in undoped 
graphene rely on the conventional phonon mediated BCS 
mechanism |14{ . which leads to an s— wave superconduc- 
tivity with low Tc values for realistic carrier densities due 
to the vanishing density of states of relativistic parti- 
cles. However, there is an alternative route to supercon- 
ductivity, wherein repulsive microscopic interactions give 



rise to attraction in a d-wave channel [l5| . This alterna- 
tive route becomes viable when graphene is doped to the 
M point of the Brillouin zone corresponding to 3/8 or 
5/8 filling of the tt band (pristine graphene corresponds 
to 1/2 filling). At this filling factor, a logarithmic Van 
Hove singularity originates from three inequivalent saddle 
points, and the FS also displays a high degree of nesting, 
forming a perfect hexagon when third and higher neigh- 
bor hopping effects are neglected [71, ll2| (FigH]). The 
combination of a singular DOS and a near-nested FS 
strongly enhances the effect of interactions [18|-|20[ , allow- 
ing non-trivial phases to emerge at relatively high tem- 
peratures, even if interactions are weak compared to the 
fermionic bandwidth W. Relevant doping levels were re- 
cently achieved experimentally using calcium and potas- 
sium dopants [l6|. Also, the new technique 17| which 
employs ionic liquids as gate dielectrics allows high lev- 
els of doping to be reached without introducing chemical 
disorder. 

Competing orders: In systems with near- nested FS, su- 
perconductivity (SC) has to compete with charge density 
wave (CDW) and spin density wave (SDW) instabilities 



sidered a charge 'Pomeranchuk' instability to a metallic 



phase breaking lattice rotation symmetry, and |25l . [26 1 



21 



At the first glance, it may seem that a system with 
repulsive interactions should develop a density-wave or- 
der rather than become a superconductor. However, to 
analyze this properly, one needs to know the susceptibil- 
ities to the various orders at a relatively small energy, 
Eq, at which the order actually develops. The couplings 
at Eq generally differ from their bare values because of 
renormalizations by fermions with energies between Eq 
and W. At weak coupling, these renormalizations are 
well captured by the renormalization group (RG) tech- 
nique. 

Interacting fermions with a nested FS and logarith- 
mically divergent DOS have previously been studied on 
the square lattice using the RG methods [18|-l21[, where 
spin fluctuations were argued to stimulate superconciduc- 
tivity. However, analysis also revealed near degeneracy 
between SC and SDW orders. The competition between 
these orders is decided by subtle interplay between de- 
viations from perfect nesting, which favor SC, and sub- 
leading terms in the RG flow, which favor SDW. In con- 
trast, the RG procedure on the honeycomb lattice un- 
ambiguously selects SC at leading order, allowing us to 
safely neglect subleading terms. The difference arises be- 
cause the honeycomb lattice contains three saddle points, 
whereas the square lattice has only two, and the extra 
saddle point tips the delicate balance seen on the square 
lattice between magnetism and SC decisively in favor of 
superconductivity. A similar tipping of a balance be- 
tween SC and SDW in favor of SC has been found in RG 



studies of some Fe-pnictide superconductors |22|.|23|. 

In previous works on graphene at the M point, vari- 
ous instabilities were analyzed using the rando m-p hase 
approximation (RPA) and mean field theory. In Il2l . the 
instability to d-wavc SC was studied, whereas [2J] con- 



considered a spin density wave (SDW) instability to an 
insulating phase. Within the framework of mean field 
theory, utilized in the above works, all of these phases are 
legitimate potential instabilities of the system. However, 
clearly graphene at the M point cannot be simultane- 
ously superconducting, metallic and insulating. The RG 
analysis treats all competing orders on an equal footing, 
and predicts that the dominant weak coupling instability 
is to superconductivity, for any choice of repulsive inter- 
actions, even for perfect nesting. Further, the Ginzburg- 
Landau theory constructed near the RG fixed point fa- 
vors the d + id state. 

The model: We follow the procedure developed for the 
square lattice [2l| and construct a patch RG that consid- 
ers only fermions near three saddle points, which dom- 
inate the DOS. There are four distinct interactions in 
the low energy theory, involving two-particle scattering 
between different patches, as shown in Figl2] 

The system is described by the low energy theory 

3 1 

^ "" X! "^^a (5r - Ek + /i) V'a - ■^fjii^ii'i'ipai'a (1) 

where summation is over patch labels a, /3 = Mi, M2, M3. 
Here £k is the tight binding dispersion, expanded up to 
quadratic terms about each saddle point. For example. 



near point Af 1 , the tight-binding model |27| predicts dis- 
persion Ek = 2'!T^aH{{Sk^)^ - VSSkJky + 0{{dk)'')), 
where t is the nearest neighbor hopping, and a is the 
lattice constant, and (5k ~ k — k^/i . The chemical poten- 
tial value fi = describes system doped exactly to the 
saddle point. We note that while the existence of saddle 
points is a topological property of the FS and is robust 
to arbitrarily long range hopping, the FS nesting is spoilt 
by third and higher neighbor hopping effects [7|, |l2| . In- 
equivalent saddle points are connected by a nesting vec- 
tor Qci/3 = Cq — e^ (FigH)). A spin sum is implicit in the 
above expression, and the interactions are assumed to 
be spin independent. The short-range interaction model, 
used in our analysis, is expected to provide a good ap- 
proximation under the conditions of metallic screening 
arising due to the states near the FS. We further assume 
that screening is insensitive to the level of doping rela- 
tive to the M point. While these assumptions introduce 
a large uncertainty into the bare values for the interac- 
tions, we will show that precise knowledge of these bare 
values is not required for determining the final state. 

The patch structure of the interactions is restricted by 
momentum conservation, which allows only the four in- 
teractions in ([T]). The Umklapp interaction (73 is allowed, 
because it conserves momentum modulo a reciprocal lat- 
tice vector. All four interactions in ([T} are marginal at 
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FIG. 2: Possible interactions in the patch niodeL (a) Feyn- 
man diagrams representing aUowed two-particle scattering 
processes among different patches, Eq[T] Solid and dashed 
lines represent fermions on different patches, whereas wavy 
lines represent interactions, (b) Pictorial representation of 
these scattering processes, superimposed on a contour plot of 
the energy dispersion. Each scattering process comes in three 
flavors, according to the patches involved. However, it follows 
by symmetry that the scattering amplitudes are independent 
of the patches involved, and therefore we suppress the flavor 
labels. 



tree level, but acquire logarithmic corrections in pertur- 
bation theory, which come from energy scales i? < A, 
where A w t is the energy scale at which higher order 
corrections to the dispersion become important. 

Logarithmic divergences in perturbation theory analy- 
sis indicate that the problem is well suited to study using 
RG. The building blocks of the RG are the susceptibilities 
in the particle-particle and particle-hole channels, IIpp 
and Ilph, evaluated respectively at momentum transfer 
zero and at momentum transfer Qa^p (FigH)). Similarly 
to [2l[, we have 
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npp(0) = ^ln 



max (T, ^) T ' 



(2) 



Hph Qq//3 = -r In 777^^ In 777^ TT: 

4 max(J,^j max(J ,/.t,i3) 

and Ilph (0), IIpp (Qq^/j) = ^0 In i„a^|-r^^) , where A is our 
UV cutoff (FiglD) and T is the temperature. The single 
spin density of states at a saddle point is i/q In ^ax (t u) ' 
The additional log factor in npp(O) (Cooper channel) 
arises because ek = f-k, generic for any system with 
time reversal or inversion symmetry. In contrast, the ad- 
ditional log factor in nph(Qa5^^) arises from nesting of 
the FS, and is cut in the IR by any term that spoils the 
nesting, such as third neighbor hopping ^3 or doping /i 
[12l |. We assume max(i3,//) ^ A, so nph(QQ^;3) and 
npp(O) are of the same order under RG. 

RG equations: The RG equations are obtained by ex- 
tending the approach developed for the square lattice 



problem [20[ to the number of patches n > 2. The num- 
ber of patches matters only in diagrams with zero net 
momentum in fermion loops, since it is only there that 
we get summation over fermion flavors inside the loop. 
The only zero-momentum loop with a log divergence is 
in the Cooper channel. Moreover, only the (73 interaction 
changes the patch label of a Cooper pair, therefore, the 
number of patches affects only diagrams where two 173 in- 
teractions are combined in the Cooper channel. With log- 



arithmic accuracy, using y = npp(k 



as the RG time, we obtain the /3 functions 
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(3) 



Here ^1(2;) = dnpM)/du « nph(Q)/npp(0) is the 
'nesting parameter' J2Ct |2l|. This quantity equals one 
in the perfectly nested limit. For non-perfect nest- 
ing, di{y) has the asymptotic forms di{y = 0) = 1, 
di{y :^ 1) = ln|A/i3|/^, and interpolates smoothly 
in between. Since the RG equations flow to strong cou- 
pling at a finite scale j/c, we treat < di(i/c) < 1 as a 
parameter in our analysis. 

The /3-functions, Eq.Q, reproduce the two-patch RG 
from [20| when we take n = 2, and negl ect subleading 
O(log) divergent terms ((^2, 3(2/) from [20|), and also re- 
produce for n — 2 the RG equations for the Fe-pnictides 
[231 . Graphene near the Van Hove singularity however is 
described by n = 3. 

We note from inspection of ([3]) that 51 , 52 and 53 must 
stay positive (repulsive) if they start out positive. This 
follows because the (3 function for 52 is positive definite, 
and the /3 functions for 51 and 53 vanish as the respec- 
tive couplings go to zero. However, 54 decreases under 
RG and eventually changes sign and becomes negative. 
As we will see, 53 — 54 becomes large and positive under 
RG, driving an instability to a superconducting phase. 
However, the positive 53 coupling penalizes s-wave su- 
perconductivity, so pairing occurs in a higher angular 
momentum (d-wave) channel. 

We integrate our RG equations with n ~ 3 from start- 
ing from 5i = 5o = 0.1 and modeling di as di{y) = 
l/-\/l -I- y. The results are plotted in Fig 13] Similar re- 
sults are obtained if we just treat di as a constant. The 
couplings diverge at a scale j/c ~ l/SOi corresponding to 
a critical temperature and ordering energy scale 



Tc, £^0 ~ Aexp(-A/^5oi/o)- 



(4) 



Here A is a non-universal number that depends on how 
we model di{y). For di = 1 (perfect nesting, correspond- 
ing to zero third neighbor hopping ^3), we obtain A = 1.5. 



An RPA-typc estimate of go is outlined in the online sup- 
plementary material. While Tc and Eq are exponentially 
sensitive to go, thus introducing a considerable uncer- 
tainty to our estimate, a strong enhancement of charac- 
teristic energy scales relative to the BCS result is evident 
from Eq.©. 

A similar y^ dependence arises in the treatment of 
color superconductivity 28| and in the analysis of the 
pairing near quantum-critical points in 3D [29i |. It re- 
sults in a Tc that is strongly enhanced compared to the 
standard BCS result Tc ~ exp{—A'/goi'o)- It should be 
noted that the enhancement of Tc in (j4]) arises from weak 
coupling physics. It is distinct from the high Tc supercon- 
ductivity that could arise if the microscopic interactions 
were strong Refs. 13l.l30l432|. 

Returning to our RG analysis, we note that near the 
instability threshold, 51,52,53 — ^ 00 and 54 — )■ —00, with 
—54 > 53 > 52 > 5i- This observation may be made 



as 
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G, 



Vc-y 



(5) 



Substituting into Eq|31 we obtain a set of polynomial 
equations, which may be solved for the co-efficients Gt 
as a function of di{yc). The solution is plotted in the 
inset of FigH Note that -G4 > G3 > G2 > Gi for 
all values of di{yc) satisfying < di{yc) < 1. We have 
verified that any choice of repulsive bare couplings leads 
to the same limiting trajectory (see online supplementary 
material) . 

Susceptibilities: We now investigate the instabilities of 
the system by evaluating the susceptibilities x for various 
types of order. To analyze the superconducting instabil- 
ity, we introduce infinitesimal test vertices corresponding 
to particle-particle pairing into the action, C = Co + SC, 
where £0 is given by ([T]) and 
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RG scale y 



FIG. 3: Flow of couplings with RG scale y, starting from re- 
pulsive interactions. Note that the coupling p4 changes sign 
and becomes attractive, leading to a (superconducting) in- 
stability at the energy scale yc (Eq[4|). Inset: Critical cou- 
plings Gi (O near yc as a function of the nesting parameter 
at the ordering energy scale, d\(yc). The dominance of su- 
perconductivity over spin density wave order arises because 
—Ga > G2 for all values of di(yc). Initial conditions: The RG 
flow is obtained by numerical integration of (O with initial 
conditions gi{0) =0.1, and modeling the nesting parameter 
as di(y) = 1/^1 -I- y. The qualitative features of the flow 
are insensitive to initial conditions, and to how we model d\. 
The critical couplings (inset) are universal, and independent 
of initial conditions. 



precise by noting that close to yc, the interactions scale 



a=l 



one test vertex for each patch. The renormalisation of 
the test vertices is governed by the equation [20[ 



(7) 



which can be diagonalized by transforming to the eigen- 
vector basis 
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Here Ac is an s-wave order, whereas Aq and A;, corre- 
spond to order parameters that vary around the Fermi 
surface as Acos(2(/3) and Asin(2iy9), where tp is the an- 
gle to the X axis (see Fig |4]). Such dependence de- 
scribes d-wave superconducting orders (SCd), since the 
gap changes sign four times along the FS. In 2D notation, 
the two order parameters A^ and A^ correspond to dxy 
and dx^-y^ superconducting orders respectively. 

Notably, we find the s-wave vertex Ac, Eq.®, has a 
negative eigenvalue and is suppressed under RG flow ^ . 
This is to be expected given that we started out with 
repulsive microscopic interactions. At the same time, the 
d-wave orders A^ and A;, have (identical) eigenvalue 53 — 
54, which may be negative at the bare level but definitely 
becomes positive under RG, indicating an instability in 
the d-wave channel. We solve ([7]) for the d-wave orders 
by substituting the scaling form of the interactions ([5]) , 
and find that the d-wave susceptibility diverges near yc 



as 

xsc,{y)^^^-{y.-yf^^^-^^\ (10) 

where, we remind, G^ — G^ > 0. 

The divergence of the SCd susceptibihty indicates an 
instabihty to d wave superconductivity under RG, with 
the Kcos{2if) and Asin(2(^) order parameters having 
identical susceptibihty. The degeneracy of the two d- 
wave orders is guaranteed, since the dx2_y-2 and dxy func- 
tions belong to the same two dimensional irreducible rep- 
resentation of the lattice point group 12,ll3|. However, 
this does not guarantee that d-wave superconductivity 
will develop, since the SCd instability must compete 
against the tendency for density wave formation. 

To investigate density wave formation, we introduce 
test vertices representing pairing of particles with holes 
on a different patch. The particles and holes may pair in 
the charge channel, forming CDW, or in the spin chan- 
nel, forming SDW. We compute the rcnormalization of 
the pairing vertices under RG, and find that the CDW 
vertex is suppressed by interactions, but the SDW vertex 
^SDW is enhanced, similar to [20] ■ The SDW suscepti- 
bility XSDW diverges near j/c as 
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This describes a potential instability towards SDW for- 
mation, which will compete with the SCd instability. The 
SDW instability arises provided there is at least partial 
nesting i.e. the nesting parameter di{yc) ^ 0. However, 
since -G4 > G2 for all < di{yc) < I (FigE] inset), it 
follows from comparison of Eg ITT] and EqlTU]that the SCd 
susceptibility diverges faster than the SDW susceptibil- 
ity, for all values of nesting. At perfect nesting {di = 1), 
the SCd susceptibility diverges as (j/c — y)~^'^, whereas 
the SDW susceptibility diverges only as {yc — y)~^- As we 
move away from perfect nesting, the SCd susceptibility 
diverges faster, and the SDW susceptibility diverges more 
slowly, so that SCd is the leading instability for all values 
of nesting, within validity of the RG. This is in contrast 
to the square lattice [20[, where at perfect nesting the 
SDW and SCd instabilities have the same exponent un- 
der RG, with subleading terms lifting the degeneracy in 
favor of SDW, which is in turn overtaken by SCd at some 
di <1. 

We also considered the possibility of ordering in a chan- 
nel exhibiting only a log divergence e.g. the Pomer- 
anchuk ordering. However, we found that such orders 
cannot compete with superconductivity (see online sup- 
plementary material). Finally, the phonon-mediated at- 
traction in the pairing channel could induce s-wave su- 
perconductivity provided that it overwhelms the elec- 
tronic repulsion in the s-wave channel at the Debyc fre- 
quency scale, loj:) < A. However, the s-wave coupling 



(2g(3) -f- g(4)) remains positive and grows ever more re- 
pulsive under our log RG. Thus, as long as the nesting 
parameter is not too small (2G3 > — G4 for di > 0.05), 
the s-wave pairing appears to be unlikely to win. 

Competition of d-wave orders below Tc : We now inves- 
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FIG. 4: Possible superconducting orders that could develop 
at the M point, (a) A d^2_y2 or d^^y state would be realised 
if K2 < in the Landau expression for the free energy, Ea ll2l 
(b) The d^2_y2 and d^y orders can co-exist if _R'2 > in Eq I12I 
A microscopic calculation indicates that the states (b) have 
lower free energy. 

tigate the competition of the dx2-y2 and dxy supercon- 
ducting orders ([8]) below Tc- In this regime, the system 
may cither develop one of these two orders, or a linear 
combination of the two. The ordered state that mini- 
mizes the free energy wins. The hexagonal lattice point 
group symmetry dictates that the free energy below Tc 
must take the form ISSl 



F = a{T-Tc)i\Aa\' 
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-K,\Al + Al\' + 0{A'^) 
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with Ki > 0. This free energy allows for two possi- 
ble superconducting phases. If K2 < then a dx2-y2 
or a dxy superconducting state would arise, whereas if 
if K2 > then the dx2-y2 and d^y orders can co-exist 
[33| . We now calculate K2 microscopically (an alterna- 
tive but equivalent microscopic treatment is provided in 
the online supplementary material). 

We begin by writing the free energy as the sum of the 
free energy on three patches, 

F = F{A,) + F{A2) + F{As), (13) 

where the free energy on a patch is given by the standard 



Landau expansion 

F(A0 = a'(T-Te)|A,|2 + /^|A/, A' > (14) 

In this expression, it is essential to realize that Ai, A2, 
and A3 are not independent, but must be expressed in 
terms of the two parameters Aq and A;,, Eq.®. Rewrit- 
ing ()13p and (J14p in terms of the two independent vari- 
ables Aa,b, we obtain Eq. (fT2| with Ki = ^K > and 
K2 = hK > 0. This implies the co-existence of d^2_y2 
and dxy orders. Minimization of the free energy (J12p with 
K2 > leads to |Aa| = |Af,| and Arg(Aa/A{,) = 7r/2. 
This order parameter can be rewritten as a three compo- 
nent vector in the patch basis, which takes the form 



[r 



A, ± iAfc = A(l, e±2W3^ gT2W3-)_ 



(15) 



This corresponds to a superconducting gap that varies 
around the FS as A exp(±2i(^). Such an order parameter 
corresponds to d+id (or d—id) superconductivity (FigUl), 
and is a spin singlet analog of the p + ip state that has 
been predicted for Sr2Ru04. 

Conclusions: The robustness oi d + id superconductiv- 
ity in the weak coupling limit, demonstrated by the above 
analysis, leads us to believe that the graphene based 
chiral superconductivity can be realized experimentally. 
While our analysis is controlled for the weak short-range 
interaction model, several questions pertaining to the be- 
havior of realistic systems should be clarified by future 
work. Determination of the phase structure for inter- 
actions of moderate strength and of long-range character 
remains an open problem, as docs an accurate estimate of 
Tc and the role of disorder, against which d-wave super- 
conductivity is not protected. The graphene based d-\-id 
superconductivity, if realized in experiment, will play a 
vital role in the development of technology designed to 
exploit topological superconductivity. 
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Now, we start with the Coulomb interaction V{q) = 
27re^/(ga^) (we keep track of the lattice scale, so that 
the Coulomb interaction in momentum space continues 
to have dimensions of energy). Let us integrate out all 
states between the bandwidth W and the UV cutoff A. 
We will take into account the effect of these states in 
the RPA. The RPA can be formally justified by appeal- 
ing to the large number of fermion flavors present in the 
problem, N=6 (3 patches and 2 spins) . The effect of the 
high energy states in the RPA is to screen the Coulomb 
interaction 



V{q,uj) 



1 



27re 



'f;^NU^iq,u;) 



(18) 



where A^ = 6 is the number of fermion flavors participat- 
ing in screening and IIa is the single species polarisation 
function obtained by integrating out states with energies 
greater than A. The interactions gi and 34 do not alter 
the patch label of electrons - treat these as having mo- 
mentum transfer q = qa ~ l/(7ra), where tq\a?' « A. 
Meanwhile, the interactions 52 and gz transfer momen- 
tum Q = 1/a, equal to the nesting vector connecting 
in-equivalent patches. Neglecting any frequency depen- 
dence of the interaction, 



5i(0)«.g4(0)«y(9A,0); 
ff2(0)«.93(0)«l^(Q,0) 



(19) 
(20) 



SUPPLEMENT 

Estimating Tc 

We start by noting the tight binding dispersion 

E = i^A/l + 4cos(7rfc2;a)(cos(7rfc3;a) -t- cos(v37rfcj,a)) 

(16) 
where t = ieV is the nearest neighbor hopping for 
graphene and a k, lA is the lattice constant. This de- 
fines a bandwidth W « 9ew. The linearization of the 
dispersion about the M point is valid for irka <^ 1, which 
corresponds to an energy window of width Aq ~ t. In 
this energy window, the density of states takes the form 



v{E) wi^ohiAo/^ 



(17) 



with z/Q ~ l/(27ri). 

The RG is performed starting from some UV scale A. 
Since the RG is performed with the linearized dispersion, 
it is essential that A < Aq. The natural choice would be 
to take A = Ap = i. For technical reasons, however, it 
will prove convenient to take A to be close to but slightly 
smaller than Aq. For the present, we simply treat A as a 
free parameter. 



Now, IIa should be like the polarisation function for 
an insulator with bandgap A. In particular, for large 
momentum transfer q > qA, there should be metallic type 
screening, with IIa (g > Sa) ~ ^i-^) = i'oln(Ao/A). Thus, 
we get 



gii^o = giVa = 



■ne^/{at) 



giJ^o 



9A^Q 



l + iV7rln(Ao/A)e2/(ai) 

eV(at) 
l + iVln(Ao/A)eV(ai) 



(21) 
(22) 



where we have taken vq = l/27ri. Now e^ /{at) sa 5 for 
a = lA and t = 3eV, so we are close to unitarity. In this 
limit, we have 



1 



gii^o ~ 921^0 ~ 9d,vo w giVf) w g^v^ 



iVln(Ao/A) 



(23) 



Let us take IuAq/A = 1, i.e. A = leV. There is ad- 
mittedly an arbitrariness in this choice. However, this 
arbitrariness reflects itself only in an 0(1) uncertainty in 
the prefactor for our eventual expression for T^. Calcu- 
lation of this prefactor is beyond the scope of RG, which 
can only calculate the exponent in the expression for Tc 
(logarithmic accuracy). Thus, with logarithmic accuracy, 
we have gi = go = 1/6. We take these to be the bare cou- 
plings at the UV scale A = leV, and henceforth do RG. 



Substituting into the results obtained by integrating the 
hi RG equations, we obtain an estimate 



Tc = Aexp(-1.5/VffoiA)) ~ SOOiiT 



(24) 



(up to pre- factors of order unity). If true, this would 
be a remarkable result, exceeding the critical tempera- 
ture of all other known superconductors. However, we 
have made some strong approximations in obtaining this 
result. For example, we have completely neglected sin- 
gle log terms in the RG equations. While such single 
log terms are formally subleading in the limit of weak 
coupling, they may well affect Tc. Moreover, since Tc 
is exponentially sensitive to go, our approximate calcu- 
lation of go carries an exponentially large uncertainty in 
the value of T^. 



The fixed point trajectory 

Here, we address the question of how large is the basin 
of attraction for the fixed point investigated in the main 
text. We show that the basin of attraction for the fixed 
trajectory includes the entire parameter space of weak 
repulsive interactions. We recall the RG equations 



dgi 
dy 
dgs 
dy 
dgi 

dy 



2rfi.gi(g2 -51), 



dg2 
dy 



diigl + gl), 



-{n - 2)gl - 2(7354 + 2^133(252 - 5i), (25) 
-{n- l)gl- gl 



We note that the equations ([25]) are homogenous, and the 
j3 function for 172 is positive definite. If we assume that 
the initial 172 interaction is positive (repulsive), then 172 is 
monotonically increasing under RG, and can be treated 
as a proxy for the RG time, following [3J]. Making the 
substitutions gi = Xig2, 53 — 3^352 and 174 = Xig2, we 
can rewrite (PS|) for n = 3 as 



-xi 



-X3. 



-X4 



2x1(1 — a;i) 

1 + 4 ' 

2^1X3(2 — Xi) 



2X3X4 



di{l+xl) 



, (26) 



2x3 



di{l + xl) 



dxi 
d In 52 

din 32 

dx4 
din 52 



The fixed points of ([26l) (e.g., solutions with constant xi, 
X3, and X4) correspond to fixed trajectories of the RG 
flow. When solving (PS| with dxt/dlng^ = 0, di should 
be interpreted as di{yc), and we should restrict ourselves 
to solutions with real values of xi, X3 and X4, with xi > 
and X3 > 0. The latter constraint follows because the (3 
functions for gi and 53 ([25]) vanish when the respective 
couplings go to zero, and so gi and 173 cannot become 
negative if they start out positive. 



The set of non-linear algebraic equations for Xi reduces 
to 7th order equation on. say, xi, hence in general there 
arc 7 different fixed trajectory. However, three of them 
correspond to negative values of either xi or X3, and three 
fixed trajectories are unstable, as we verified by solving 
the set (P5|) near the fixed trajectory. This leaves the 
fixed trajectory discussed in the main text as the only 
stable fixed point of ((26|) that is compatible with the 
above constraints. Thus, any choice of weak repulsive 
interactions leads to the same fixed trajectory. 

The solutions for Xi along the fixed trajectory can be 
obtained analytically if we assume that the bare value of 
the exchange coupling 171 is zero, in which case 171 = 
holds during RG fiow, and xi = 0. The set of two alge- 
braic equations for X3 and X4 at the fixed point then re- 
duces to 4th order polynomial algebraic equation, which 
can be solved exactly. Out of 4 soltions, two correspond 
to imaginary X3 and one to a negative 173. This leaves only 
one fixed trajectory, consistent with intial conditions. 



Ordering in 0(ln) divergent channels 

Here we consider the possibility of ordering in an 0(ln) 
divergent channel, and show that it cannot compete with 
superconductivity. First, we recall the scaling form of the 
superconducting susceptibility. 



Aa,fe(y) 

A,,,(0) 



= xscdiy) ^ ivc-y) 



_„^2(G4-G3) 



-go 



gi{y) 



2(G4-G3) 



(27) 



We wish to contrast this with the susceptibility in an 
0(ln) divergent channel. We therefore introduce a ver- 
tex corresponding to particle-hole pairing on the same 
patch, and examine how it renormalizcs under RG. We 
find a scaling solution for the susceptibility, generic for 
any ordering in a 0(ln) divergent channel, which takes 
the form 



X^(2/c-2/)"/^^(3o/.94(y))' 



/go 



(28) 



where a is some linear combination of the Gi with 0(1) 
coefficients. Naively, such susceptibilities will also di- 
verge as y — )- j/c if a < 0, although the exponent will be 
parametrically smaller than pT]) by ^/gH- However, we 
argue that not only is the exponent for these divergences 
parametrically small, but in fact such divergences lie out- 
side the range of justifiable applicability of the RG. To 
understand why, it is essential to remember that the one 
loop RG only applies upto an energy scale j/i when the 
couplings become of order one. (The limiting scale j/i 
may actually be even smaller once we take into account 
self energy I](w, Aif) oc g^wlog^ A/|a;| (Ref.^)). 

At the scale yi, ([28]) gives x(2/i) = exp (ay^lngo)-In 
the weak coupling limit, go — )> and xivi) ~ 1- There- 
fore, the susceptibility in a In divergent channel is not 



significantly enhanced within the region of apphcabihty 
of the one loop RG. In contrast, for In divergent chan- 
nels like SCd, 1^^ gives Xscdivi) ^ 3o ^^'', which goes to 
infinity as go ^ 0- Therefore, only susceptibilities in In 
divergent channels are strongly enhanced in the regime 
of justifiable applicability of weak coupling RG. 



Hubbard-Stratonovich treatment of 
superconductivity 



We have written the action in a Gorkov-Nambu spinor 
form, introducing the particle and hole Green functions 
G_|_ and G'„. These Green functions are diagonal in 
Fourier space, and have the form G^^ (w, k) = iuj^ (e^ — 
/i) where w is a Matsubara frequency, £k is the energy of a 
state with momentum k and /i is the chemical potential. 
We now integrate out the fermions in Eql33]to obtain an 
exact action in terms of the order parameter fields alone. 
This action C'iAa, A* , Aj,, A*) takes the form 



Here, we provide details of the Hubbard-Stratonovich 
treatment used to investigate the superconducting phase 
at temperatures lower than Tc. We begin by writing the 
partition function in the path integral formalism as Z = 

/I?['0, ?/>] exp(— j£[-0,-0]), where 
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^i'iidT-£-k + ^J■)i^a 
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53 53 54 _ 




V'3V'3 



(29) 



(30) 



Here a is a patch index and the momentum, frequency 
and spin indices have been suppressed for clarity. We 
have included only the 'pair hopping' interactions 53 and 
34 since these are the only interactions that contribute 
to d-wave superconductivity. 

As discussed above, 173 and 174 flow under the 3-patch 
RG so that 



53 - 54 = A > 0. 



(31) 



When this is the case, then the interaction matrix in 
Ea. (P5|) has two eigenvectors with degenerate negative 
eigenvalues. These reflect the two possible d-wave su- 
perconducting phases, which have identical instability 
threshold. We introduce two 3x3 matrices in patch 
space, di and ^2, where 



1 



V2 



diag(0,l,-l): 4 



-diag(l, — , 
3 **^ ' 2' 



1 



2' 
(32) 

These matrices obey Ti (di) ~ 1, Ti (d^) = 1 and 
Tt {dadb) ~ 0, where the trace goes over the patch 
space. Using these matrices, we can define the or- 
der parameters of the two d-wave instabilities as 
Aa = 2X{tpdafp) and A;, = 2A('(/'rf6^). We can 
now decouple the quartic interaction in Ea l29l us- 
ing a Hubbard Stratonovich transformation, and 
can hence rewrite the partition function as Z ^ 
J D[^P, ^, Aa^t, a;,] exp(- / C'l^-j, V-, Ai, A*, A2, A*]), 
where 



£' = 



+ 
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G-^ Aada + Abdb ' 
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i^f^l 




[ Alda + A*4 G-_^ \ 
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-\^b? 




(?, 



C = Tr In 



G- 

Alda + Aldb 



Aada + Abdb 

GZ' 



lAal 



4A 
(34) 
The trace goes over patch space and over the spinor 
space. We expand this action in small Aa^b to order 
A^, exploiting the fact that the Green functions com- 
mute with the order parameter matrices, and the trace 
over patch space vanishes for any expression with an odd 
number of da or db matrices. We make use of the iden- 
titles Tr (dl) = Tridl) - 1, Tv (d^) - Tr(4) = 1/2, 
Tt {d^dl) = Tt (dadbdadb) = 1/6, transform from parti- 
tion function to free energy and obtain, up to an overall 
factor. 



F 



IAJ2) 



4A 



Tr (G+G_; 



(35) 



-A' ( IAJ4 + lAJ^ 



\Aa\'\Ab 



AlAl + AlA 



ha 



where if = Tr {G+G-G+G-) > 0. 

Superconductivity sets in when the coefficient of the 
quadratic terms first becomes negative, which leads to 



Tc~ A'exp(-1/VA). 



(36) 



The nature of the superconducting phase below Tc is 
controlled by the anisotropic quartic term. Since K > 0, 
minimization of the quartic term leads to |Aa| = |A;,| 
and Arg(Aa/Ab) = tt/2. The full superconducting order 
parameter is thus 



Aa±^Ab = AJ-[l,--±^-,--T^- 




(l,e±2W3^eT2W3y (37) 



4A 



